Jet Quenching in High Energy Heavy Ion Collisions 
by QCD Synchrotron-like Radiation 



(N 
O 

o 



> 

m 
o 

<N 
O 

i 

Q . 
^ ■ 



X 



E.V. Shuryak and I. Zahed 
Department of Physics and Astronomy, 
State University of New York, Stony Brook NY 11794-3800, USA 
(February 1, 2008) 

We consider synchrotron-like radiation in QCD by generalizing Schwinger's treatment of quantum 
synchrotron radiation in QED to the case of a constant chromomagnetic field. We suggest a novel 
mechanism for jet quenching in heavy ion collisions, whereby high-pt partons get depleted through 
strong (classical) color fields. The latters are encountered in the color glass condensate or in the form 
of expanding shells of exploding sphalerons. Unlike bremsstrahlung radiation through multiple soft 
rescattering, synchrotron radiation converts a jet into a wide shower of soft gluons. We estimate the 
energy loss through this mechanism and suggest that it contributes significantly to the unexpectedly 
strong jet quenching observed at RHIC. 



I. INTRODUCTION 
A. Radiation in Various settings 

By synchrotron-like radiation we mean radiation emit- 
ted by a charge moving in an external field that is strong 
enough not to allow for a perturbative treatment. The 
strong external field problem requires the exact trajecto- 
ries classically, and the exact propagators quantum me- 
chanically. Throughout, we will refer to synchrotron ra- 
diation as the part of the radiation stemming from within 
the strong field region, while the radiation from the out- 
side field region (if present) will be still referred to as 
bremsstrahlung radiation. For ultrarelativistic particles 
of energy E = 7 m, bremsstrahlung radiation is usually 
collimatcd along the particle trajectory in two cones of 
opening angle I/7 <C 1, while synchrotron radiation di- 
verges away from the particle trajectory. 

In QED synchrotron radiation usually takes place 
within a magnet as illustrated in Fig. |l|a. The classi- 
cal trajectory includes the circular bending between the 
incoming and outgoing straight lines, and the radiation 
is emitted tangent to the arc length. Photons move on 
straight lines since in QED they are not affected by the 
external magnetic field. 

Another interesting case of classical synchrotron-like 
radiation is that of an ultrarelativistic rotating charge in 
a strong gravitational field such as the one encountered 
near the horizon of a black-hole (and subsequent liter- 
ature). In this case, both the charge and the photon are 
gravitationally deflected as illustrated in Fig. [j]b. The 
result is a significant reduction of the radiation loss: the 
total radiation yield is reduced by a factor of 7 2 in com- 
parison to the yield from standard synchrotron radiation 
for the same curvature. Also, the radiation length for 
each particular direction is actually the entire circle, not 
just an arc length of order I/7 as in the magnet. 

In this paper we consider an ultrarelativistic charge in 
QCD (parton, quark, gluon) going through a constant 
chromomagnetic field as illustrated in Fig. [j]c. The mo- 
tion of the initial charge and the ensuing radiation are 



both strongly affected by the chromomagnetic field. If 
classical geometrical optics can be used g] (and for hard 
jet quenching it is not the case) the classical motion of 
the external particle and the radiation is described by 
non-trivial trajectories. For soft radiation, the emitted 
gluons are bent along circular paths of smaller curvature 
radii. As a result, part of the radiation gets trapped 
in the near field region and never makes it to infin- 
ity. This makes standard radiation calculations obsolete. 
Moreover, the radiated gluons carry different charges 
and therefore move in different directions. The QCD 
chromomagnetic field resolves both momentum and color 
thereby acting as a double (squared) Newtonian prism. 
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FIG. 1. Schematic representation of synchrotron- like ra- 
diation in three cases: (a) in a magnetic field; (b) in a grav- 
itational field (rotating charge around a black hole); (c) in a 
chromomagnetic field. 

Another instructive classical description to keep in 
mind is the instantaneous distribution of the radiation 
field in the ultrarelativistic case. In QED the nonzero 
field strength is mostly located around a spiral-like curve 
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known as the evolute of the circle. Since both the charge 
and the photon moves light-like, the distance from the ra- 
diation point to the charge is the same as to the photon 
along a straight line. Details of the field distribution for 
QED synchrotron radiation can be found in || . Its ana- 
log for QCD in the simplest abelian-type external field 
when the radiated gluons undergo planar rotation has a 
cycloidal shape. The radiation field is the same at any 
instant of time but moves (together with the charge) at 
the speed of light. 

B. Synchrotron Radiation versus 
Bremsstrahlung Radiation in QED 

The number of soft prompt gluons per frequency lu 
expected from bremsstrahlung and synchrotron radiation 
can be estimated classically. The Lienard expression for 
the power radiated from a dipole 

3m 2( ds ' ( ) 

was derived in 1898. It holds even ultrarelativistically 
since the power (energy over time) is relativistically in- 
variant. The acceleration in ([!]) relates to the field 
strength F and the particle energy E by the Lorentz for- 
mula. Classically, the energy loss depends quadratically 
on both E and F, i.e. P ~ e 4 F 2 E 2 /m 4 . This result does 
not hold at very large energies though. Indeed, it only 
holds in the range m 

e 3 (F/m 2 )(E/m) < 1 , (2) 

where we have separated 3 dimensionless physical fac- 
tors: the coupling constant, the field strength in units of 
m 2 and the relativistic gamma factor. For classical QED 
synchrotron radiation in accelerators, this condition of 
course holds even with 7 = E/m ~ 10 5 (LEP). The rea- 
son is the small electromagnetic coupling and the small 
magnetic fields H/ra 2 <C 1. 

In the QCD problem we are interested in, the fields are 
strong, the external charge is massless and the induced 
radiation recoil is large. The condition (||) does not hold, 
and we have to reassess the power radiated from first 
principles. However, we expect that qualitatively the 
differential spectra to follow general lore at small fre- 
quencies. In particular, the number of bremsstrahlung 
photons per frequency should be of order 

dN B a 1 

— - r*j 

du> tt u> 

while the number of synchrotron photon per frequency 
(per orbit) should be of order 

cWs _ a (lu/luq) 1 / 3 
dio it u> 

where luq = eH/E is the synchrotron frequency with E 
the energy of the external charge. At very low frequencies 



bremsstrahlung of course dominates, but not at large fre- 
quencies. Integrating (|j) up to some maximal frequency 
yields 

N S ~ - (uw/w ) 1/3 (5) 

TT 

which amounts to an energy loss AE$ ~ awma S /wJ' 3 . 
The maximum frequency is bounded by the classical 
characteristic frequency lu c = 37 3 o->o- As a result, the 
maximum number of synchrotron photons emitted per 
cycle is Ns ~ cry, which is seen to grow linearly with the 
energy. The maximum synchrotron energy loss AE$ ~ 
a 7 4 u>o grows cubicly with the energy. In comparison, the 
power following from the Lienard expression ([I]) grows 
quadratically with the energy AE ~ a E 2 , and since 
one circle takes time proportional to E, both expressions 
agree. 

The above qualitative reasoning is entirely classical. In 
quantum mechanics the photon back reaction cannot be 
ignored, especially for large energy losses with lu ~ E (in 
units c — h — I). Although in modern accelerators like 
LEP 7 r~j 10 5 resulting into an enhancement in radiation 
per frequency that is about 15 orders of magnitude high, 
the photon energies are usually much smaller than the 
energy of the rotating particle and the back reaction can 
be ignored. 

In the QCD case to be considered below, this is not 
true. Using quantum cutoff to max — E <C to c , we then 
find that the number of bremsstralung photons grows 
logarithmically with the energy Nb ~ a In (E/uj m i n ) 
while the energy loss grows linearly with the energy 
AEb ~ aE. These results are in total agreement with 
the ones obtained through standard quantum calcula- 
tions such as Feynman graphs. For synchrotron radiation 
a cut at tt> max = E yields 

N s ~a(E 2 /eH) 1/3 (6) 

and the energy loss for a fixed width of the field region 
Az (rather than a circle) ^] 

AE s ~e 2 {eH) 2 ' 3 AzE 2 ' 3 , (7) 

which is seen to grow less than linearly. As a result, the 
relative energy loss by synchrotron radiation AEs/E ~ 
I/E 1 / 3 is found to decrease with energy, although with 
a small power. 

The present qualitative estimates are in agreement 
with the full quantum calculations to be described be- 
low. They show that bremsstrahlung radiation domi- 
nates over synchrotron radiation at large energies, al- 
though with a small relative power of 1/3. However, we 
recall that subsequent bremsstrahlung radiations at small 



*It is not small in comparison to the radiation length: Az > 
m/eH. 
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angle interfere with each other and are further reduced 
by the Landau-Pomeranchuk-Migdal (LPM) effect, while 
the synchrotron emitted quanta are lost once and for all 
without further suppression. 

C. Classical fields in heavy ion collisions 

Synchrotron-like radiation may take place in QCD set- 
tings whenever classical and strong non-abelian fields can 
be formed. We think that such strong semiclassical gauge 
fields with amplitudes A ~ l/g, naturally arise in the 
early stages of ultrarelativistic heavy- ion collisions. Cur- 
rently, there are two first principle QCD reasonings in 
favor of prompt and strong non-abelian gauge fields in 
heavy ion collisions. 

The first reasoning by McLerran and Venugopalan |t]] 
suggests that the prompt phase of a heavy ion collision 
is a Color Glass Condensate (CGC). Schematically, con- 
sider a small element Ax± of the transverse plane defined 
by the disk-shaped boosted heavy nuclei. At high energy, 
the hadronic structure functions increase at small parton 
momentum fraction x = p/E — * 0. This increase results 
in an increase in the density uq of color charges in the 
transverse 2-dimensional plane. The number of charges 
becomes ultimatly large Nq = uq Axj_ S> 1, with a to- 
tal charge of order Q ~ «/ Nq [[ The large initial color 
charge Q is at the origin of the classical color field. The 
smallest transverse size Ax± ~ l/Qs i s fixed by the par- 
ton saturation scale Q s reached when the gauge field am- 
plitude satisfies dA ~ gA 2 [0]. In this approach, the 
virtual classical fields are attached to the nuclear wave- 
function and become real after being stripped off by the 
heavy ion collision. 

The second reasoning suggests rather that the heavy 
ion collision converts the virtual classical vacuum fields 
to real fields. In other words the classical fields originate 
from the wavefunctional of the QCD vacuum. Yang-Mills 
fields have a rich topological structure in the vacuum, 
with large amplitude tunneling described classically by 
instantons [|| (for a review see (9)). Unlike the virtual 
fields close to zero strength which can be made real by 
perturbation theory, the virtual Q tunneling fields deep 
under the barrier can only be made real after an energy 
deposition sufficient to jump over the barrier. As we 
have argued recently p0| p~3|| using arguments in favor of 
the relevance of instantons to high energy hadron-hadron 
scattering in the semi-hard (soft pomeron) region 
and insights from electroweak theory p5[ , these fields are 



released in Minkowski space in purely chromomagnetic 
form, the relatives of the sphalerons in the electroweak 
theory. Indeed, much like their relatives they are classi- 
cally unstable and evolve (explode) into a thin shell of a 
coherent color field. 

The QCD vacuum is filled with relatively small in- 
stantons (and also antiinstantons), with an average size 
p 1/3 fm which is small in comparison to their rela- 
tive separation which is of the order of 1 fm. The QCD 
vacuum is characterized by the small dimensionless in- 
stanton diluteness parameter n — np 4 ~ 0.01. Therefore 
and immediatly upon their release the QCD sphalerons 
are expected to be in a dilute phase as well. When pro- 
duced in bunch like in a heavy ion collision, the QCD 
sphalerons evolve pretty much unshattered for a time of 
the order of 1 fm before they collide and get destroyed. 
More details for this process will be given below. 

In a recent paper by one of us UM the idea of jet 
quenching on coherent classical fields has been first dis- 
cussed. The importance of the synchrotron-like radiation 
is due to both of us and was briefly advertised in Jll| . 
The present paper elaborates further on this idea and 
presents detailed quantum calculations for jet quenching 
by synchrotron-like radiation. 

II. QCD SYNCHROTRON RADIATION 

In this section we proceed to estimate the QCD syn- 
chrotron radiation in a constant and Abelian-like chro- 
momagnetic field. We will derive the exact classical and 
first quantum correction in the regime uj/E < 1, and will 
provide an approximate expression for all frequencies ui. 

The problem of quantum synchrotron radiation 
in QED was addressed in a fundamental way by 
Schwinger |l6|] using the mass operator formalism. In 
this section, we extend this approach to the quantum 
synchrotron radiation in QCD. Two essential differences 
between the QED and QCD problem: i. the non-Abelian 
nature of the charge in QCD; ii. the emitted radiation 
also undergoes magnetic deflection. A quantum calcula- 
tion is required in strong chromomagnetic fields owing to 
potentially large recoil corrections, essential for large jet 
quenching. The power radiated will be sought through 
the mass operator as 

1 f dto 

--ImM QQ = / —P aa {w) (8) 

after pertinent kinematical identifications. 



'We are assuming that all partons originate from different 
nucleons and add randomly. 

*In the so called Landau method used e.g. in jl2| , those fields 
even have additional singularity, essential for the evaluation 
of the tunneling cross section. 



A. Abelian Chromomagnetic Field 

For simplicity, we consider QCD synchrotron radiation 
in a constant and homogeneous chromomagnetic field 

G^(x)=S a8 G^ (9) 
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where the abelian field strength corresponds to a constant 
magnetic field in the 3-direction, G\i = —G21 = H. The 
background gauge field associated to (0) is 



A a (x) = 6 a *A fi (x) =5 a *5^Hx 



(10) 



With our choice of the chromomagnetic background 
along the 8th color direction, the quarks and gluons can 
be diagonalized. The diagonal quarks in the fundamental 
representation carry color (a — 1,2,3) 



e a =g(T 8 
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(1,1,-2) 



(11) 



and the diagonal gluons in the ajoint representation carry 
color 



9 A = (-1) 



A 3^3 



(12) 



for A = 4,5,6,7 and g A = for A = 1,2,3,8. These 
two cases, as will be shown below, lead to qualitatively 
different radiation. The second case is basically QED- 
like. 

Quantum synchrotron radiation will be sought for 
quarks and gluons interacting to all orders in H but to 
leading order in a = g 2 /4ir between the quantized fields. 
We now present briefly the spin- 1/2 case and discuss ex- 
tensively the spin-0 case. In the semiclassical limit, both 
spins radiate at the same rate. 



B. Spin 0, 1/2 Jets 

Following Schwinger p6[ , to lowest order in perturba- 
tion theory the quark mass operator in the chromomag- 
netic field reads 



00 — iO 



ds 



M aa =ig (T ) ab (T , oa . . 

J (2tt) 4 J _ i0 cos{g A Hs) 

x e -is(k 2 -k 2 L (t&n{g A Hs)/(g A Hs)-l)) ^g A sG^ $(g A ) 

x 7^ (7 • (n 6 - k) - m)- V • (13) 
We have defined the quark 4-momentum operator as 

n 6(U (x) = id ll -e b A fl (x) (14) 
and the Bohm-Aharanov line 

:9 A H 



${g A ; x, y) = e 1 ^ (*i+wi)(*a-jBi) 



(15) 



The Bohm-Aharanov phase enforces gauge-invariance in 
the mass operator, but does not contribute to the ra- 
diation. Indeed, for an initial color-a quark emitting a 
color-b quark plus a color-A gluon, 



$(e a ;x,y) = $(e b ;x,y) x $(g A ;x,y) 



(16) 



showing that the Bohm-Aharanov line in (^3|) on the 
gluon, can be redistributed to compensate the analogue 



ones on the quarks. This procedure will be assumed 
throughout, and thereby the gluon $ contribution reshuf- 
fled. 

The occurrence of the synchrotron poles in the gluon 
propagator ( |l3| ) implies that the s- integration is infinites- 
imally shifted below the real axis in the complex s-plane. 
The prescription follows the causal prescription for the 
free propagator, 



1 



k 2 + iO 



ds < 



-is k 



(17) 



Also, since G is an antisymmetric matrix, its eigenvalues 
±iH are complex. The color precession factor is 



^g A sG 



,±is(2g A H) 



(18) 

There is a subtlety due to the positive sign in (|lj) for 
k = 0, which is the analogue of the tachyonic mode of 
a spin-1 coupled to a constant chromomagnetic field in 
the first quantized approach. This mode is at the ori- 
gin of the well-known Savvidy instability in QCD Jl7[ . 
What it says, is that in QCD the chromomagnetic fields 
themselves are in general unstable against gluon emis- 
sion. Although this phenomenon is interesting by itself, 
it clearly has nothing to do with jet ener gy l osses. 

The technique developed by Schwinger |16| can now be 
applied to (|l^) to derive the power radiated in a QCD 
synchrotron process whereby an energetic quark radiates 
through a chromomagnetic field. To avoid the unneces- 
sary algebra triggered by the spin content of the quark, 
we present the results for the spin-0 case instead. 

For a scalar quark in the fundamental representation, 
the analogue of ( |l3|) is 



M„ 



ig 2 (T A ) ab (T A ) ba 



d A k 



ds 



cos(g A Hs) 



x e -is(k 2 -k 2 L (ta,n(g A Hs)/(g A Hs)-1)) ^g A sG^ ^ <&{g A ) 

x (Ha - U b y ((H - k) 2 - m 2 )- 1 (u a - Jib)" , 

(19) 

modulo counter-terms. The arrows on FI's indicate the 
direction of the derivative. On mass-shell we expect II a ~ 
lib + k, this will hold in the classical limit. 



C. Power Radiated 

For spin-0, the power radiated follows from (||). Fol- 
lowing Schwinger |16| we obtain the chromomagnetic syn- 
chrotron emission by a scalar quark in the classical limit 
in the following form 



PaaH 



-— {T A ) ab (T A ) ba 
ir 



xlu Im 



oo-iO dr e -i(E^ A ) 2 r 3 /(24w) 

cos (Elu a t /(2ui) 



0-i0 



,m 1 2 o. _ 
X( £2 + 2 Wfc )e 



(20) 
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where in ( ]20| ) the H — subtraction is not explicitly 
shown but implied. The quark synchrotron and gluon 
rescaled frequencies are uj a = e a H/ E and uja = gA H/ E 
respectively. 

In carrying out (|2(]) the emitted gluon recoil effect on 
the jet was ignored, and so the result is entirely classical. 
We have checked that the gluon recoil effect amounts in 
the first order to the shift 



f 

U3 



(21) 



in the combination P aa /uj (the gluon multiplicity), 
thereby generalizing Schwinger's first quantum correction 
in QED to the QCD case. Of course, this substitution is 
not the complete quantum answer, but will be discussed 
below as an approximation. 

Before we discuss the complete results, let us comment 
the integrand of (^0|) . The last exponent is due to charge 
curving, and is the same as in QED. It provides a rapidly 
oscillating phase at large uj and a corresponding cutoff. 
The first exponent in the integrand of (20|) is new. It 
stems from the gluon rotation in the chromomagnetic 
field. The phase follows from the transverse contribution 
of the gluon propagator as is evident from (|l9|). The 
cosine in the denominator exhibits poles for 



Ecu a 



(2n+l) 



(22) 



which are the gluon synchrotron orbits (classically the 
gluon spin and the tachyon problem drop). The second 
contribution in (EG) is the quark synchrotron contribu- 
tion as in QED ||P 

Rewriting I/cos A as a geometrical sum of all powers 
of e lA , we may bring (pot) in the form of a sum 



uj-'P aa (uj) = -(T A ) ab (T A ) ba 

TT 



2m 2 
\[ZE 2 



x E e 

n=0 



-27T n (1— iO) 



(^(£„)+2k„K 2/3 (£„)) (23) 



with 



(! 

K n + 1 



and 



dtK 5/3 (t) , 

(l + (2n + l)(2/3)/(3£ 2 / 3 )) 3 

(1 + A 2 ) 
l + (2n+l)(2/3)/(3e 2 / 3 ) 
(1 + A 2 ) 

1/3 



(24) 






/ 3uj b \ 


X 


\2uj) 


A = 


Euj a 




LUUJ b 




2uj /TO 


3^ \~E 



(25) 



Again the quantum corrections follow from (|23|) through 
the substitution (^|]) on the RHS. The K's are modified 
Bessel functions. The sum over n in ( p3| ) sums over syn- 
chrotron orbits of width — iO except for the lowest orbit 
which is zero. It is reminiscent of the sum over 'Landau 
levels' in the Schroedinger formulation. 



D. Small and Large uj 

The preceding results are easily analyzed for large and 
small frequencies u). Since the abelian analysis with 
gA = is known from QED, we focus on the non-abelian 
part with gA =/= 0. We will show that the non-abelian 
contribution to the synchrotron radiation is strongly sup- 
pressed at small uj due to strong 'incoherence' effects pro- 
duced by the deflected radiation. At large uj the non- 
abelian contribution is equal to the abelian contribution. 
Specifically, the gA ^ contribution for small uj reads 



\9vW 

OO 

(-1)" (2n + 1) 9/4 e -|(2«+i) 3/2 V^7A 



(26) 



n=0 



which is characterized by an essential singularity at uj 
0. The gA 7^ contribution for large uj reads 



P aa (uj) « - (T A ) ab (T A ) ba 

7T 



5/2 



*(¥)(sy ^e—.-— . (27) 

At small uj the non-abelian radiation is of order 
e _ #/ a ' J 'y/cJ which is much smaller than the abelian radi- 
ation of order uj 1 / 3 . The reason is that the emitted non- 
abelian gluon brings about its own phase which strongly 
adds to the phase incoherence at small uj as is clearly 
seen in (13) and (|l9|). At large uj both the abelian 
and non-abelian radiations are comparable and of order 
w 5/2 e -2aV(3^f, 7 3 ). Indeed, we note that @ is indepen- 
dent of the gluon charge gA- For ultrarelativistic jets, the 
radiation frequencies are in the range < uj < E. Thus 
u/uj b < 7 which is way below the maximum of 7 3 . In 
light of the present observations we conclude that most 
of the jet radiation is emitted through the 'abelian' part 
of the gluon charge and in the small frequency range way 
below the synchrotron maximum since E / uj b <t^. 7 3 . The 
'non-abelian' part only enters as a correction. 



E. Energy Loss 

The 'abelian' part of the jet radiation follows from the 
summation over A = 1,2,3,8 for which gA — 0, thus 
e a = e b . In this case ( p3| ) simplifies to the QED- like 
answer 
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- (T ) ab (T ) ba 



HO (28) 



since £ n = £ and K n = 0. In this case, the total power 
emitted follows by integrating (^8|) over the gluon fre- 
quency, explicitly including the recoil effects (p3|), 



— (T )aa (T )aa r- 
7T Vo 



3~ (£corr ) 



E E 

with a different (quantum corrected) £ c 



11?/ T 



-uj/E 



(29) 



(30) 



In Fig. U we compare two spectra for some particular 
selection of parameters. The general result for (^9|) can 
be obtained by expanding around the classical result with 
the first quantum correction included 



{T ) aa (T A ) aa C (u a E*) 



2^2/3 



with the constant 

_ -tt(3/2) 8 / 3 2 5 / 3 



10sin(7r/3) T(-2/3) 



(1 



21 



0.52 . 



(31) 



(32) 



The contributions in (1 — 8/21) are the classical contribu- 
tion and the first quantum recoil correction respectively. 
So the recoil of the emitted gluon decreases the net radi- 
ation by a factor 0.62. The ratio of the contributions of 
the classical to recoil corrected spectra is 0.48. Within 
few percents these results agree with direct numerical 
estimates of the integrals including the approximate all 
orders quantum recoil corrections. This justifies a poste- 
riori the use of the first order correction in our calcula- 
tions. 

In order to further compare losses related with emis- 
sion of QED-like and QCD-like gluons, we show in Fig. || 
the leading (n=0) integral 



I 



dr m 2 
o 

[ 24w 



Euj AT 



2^3 



2E 2 



24 



(33) 



H = 1 



with the H — subtraction implied, for E = \ 
and 4 different sets of charges (sec the figure caption). 
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FIG. 2. The integral in ( |33| ) versus the radiated en- 
ergy fraction z — u/E for classical (upper figure) and quan- 
tum corrected (lower figure) cases. Four examples are shown 
corresponding to all possible charge combinations, namely 
e a = l,e b = l,g A = 3 (boxes), e a = l,e b = -2,qa = 3 
(crosses), e a = 1, e-b = = (diamonds) and 

e a =l,eb = —2,gA = (circles). 



III. APPLICATIONS TO HEAVY ION 
COLLISIONS 



In this section we will first recall for completeness some 
qualitative arguments regarding jet quenching in heavy 
ion collisions, most of which have already been discussed 



G 



in jry . We will then use the above synchrotron radiation 
results to assess jet energy loss in the color glass conden- 
sate approach and the exploding sphaleron approach. 



A. Introduction to Jet Quenching 

Jet quenching is a sort of "tomography" of the prompt 
excited system, created in high energy heavy ion colli- 
sions. Even very hard jets radiate and lose some energy 
during their passage through the system, thereby provid- 
ing information about the early stages of the collision. 

The so called quenching factor Q(pt) is defined as the 
observed number of jets normalized to the expected num- 
ber of jets calculated in the parton model without account 
for final state interactions It is usually assumed that 
hard QCD probes are under good theoretical control, and 
that one can assess the initial jet production reliably. 

Experimentally, jet reconstruction in a heavy ion en- 
vironment is very difficult to achieve. Therefore, all cur- 
rently reported results for jet quenching refer to the ob- 
served/expected ratio of the yields of single hadrons. Fur- 
thermore and in so far, large pt means pt =2-6 GeV. 
Which part of this transverse momentum comes from 
genuine high energy jets is anybody guess. The first di- 
rect evidences for jets have been recently reported by the 
STAR collaboration ]i"8| ], whereby a second particle cor- 
related with the trigger with p t > 4 GeV was observed, 
well inside a relatively narrow cone typical of jets. 

In the early theoretical studies on the subject [|l9|| , 
a rather modest jet re-scattering in the Quark-Gluon 
Plasma (QGP) has been considered. Accounting for the 
radiation effect |2(J has significantly increased expecta- 
tions for the magnitude of the result, while accounting for 
the Landau-Pomeranchuck-Migdal (LPM) effect has 
somewhat decreased the magnitude of the result. We 
will not discuss this involved subject, but only recall 
that the expected quenching factor from such studies is 
Q(p t ) =0.5-0.7 for jets with p t =10-20 GeV. 

Experimentally, a relatively modest jet quenching has 
been first observed in deep inelastic scattering for a for- 
ward jet going through cold nuclear matter (for recent 
discussion and references see [^2|). The heavy ion data 
at the CERN SPS have also shown modest quenching 
effects, but already the very first RHIC data [g3| (espe- 
cially for 7T° from the PHENIX collaboration) have shown 
that quenching of jets is very strong, with Q(pt) < 1/3. 
Subsequent discoveries that (i) at p t > 2 GeV protons 



§ Effects due to initial state interaction should be also in- 
cluded. What this means is that the parton distribution func- 
tions should be nuclear rather than hadronic, following from 
lepton-nuclei experiments. Parton rescattering in nuclei at 
the origin of the so called Cronin effect, should also be in- 
cluded in the expected yield. 



and anti-protons dominate the charge particle spectra; 
(ii) that the azimuthal asymmetry remains very strong 
even at large p t ; (in) that development of very successful 
hydro and/or cascade description of spectra without jets 
even at pt =2-3 GeV, have all led to suspect that the real 
jet quenching factor maybe even stronger. Furthermore, 
one of us even found p4 | that the very high degree of 
azimuthal asymmetry observed by the STAR experiment 
at p t =2-6 GeV [Jl8j cannot be reproduced by any amount 
of jet quenching, no matter how strong. 

One technical but important point made in the last 
paper in pl| ], is that when the quenching is strong it 
cannot be evaluated using the mean energy loss. Specif- 
ically, the quenching factor can be seen as the ratio of 
produced-and-quenched to produced spectra 



QiPt) 



J deD(e)dN/dp 2 t (p t 
dN/dp 2 (p t ) 



(34) 



where e is the energy lost in the medium and D(e) its 
normalized distribution. For small e one can expand it 
to first order, obtaining a correction proportional to the 
mean energy loss < e >= J deD(e)e. However, because 
the spectrum is so steep, this is only valid when e/pt is 
not larger than few percent. 

Therefore and for a qualitative assessment of the mag- 
nitude of the effect needed, we suggest a different simple 
approximation. Using a power parameterization of the 
spectrum 



dN 



1 



dpfipt) Pt 



(35) 



for both the observed and "hard" distributions, we obtain 
Q{p t ) = J deDieXj^—r ~ / deD(e)e~ n ^ (36) 



instead of (34). Using a simple delta- like distribution 
peaked at some fractional loss, 



D{e) = 8{e - K Pt ) 



(37) 



we have Q{pt) = 1/(1 + n) n - With n w 12 in the few- 
GeV domain at RHIC energies, a jet quenching by one 
order of magnitude would correspond to k ~ 1/4. This 
means that a mean loss of about 15-20% of the produced 
jet momentum is sufficient. However, this conclusion is 
oversimplified. As one can see from ( |36|) the quenching 
factor is dominated by small losses e/pt < 1/n ~ 1/12. 
What this means is that what matters is the probability 
to escape with as small losses as possible. 

This conclusion changes the relative role of early 
(synchrotron-like) versus late (multiple bremsstrahlung 
with LPM) effects. While the latter can be very large 
for specific geometry (LPM energy loss |2l[] is AE ~ L 2 
where L is the path in matter), the integral d36| ) would be 
dominated by surface emission with small L ~ 1 fm or so. 
Early effects emphasized in this paper, even if producing 
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less average losses, are expected to affect the probabil- 
ity D(e) at small e, preventing easy escape of some jets. 
Detailed numerical simulations (which are well beyond 
the limits of this work) are needed to understand their 
relative role. 

Summarizing this subsection, we say that the tradi- 
tional approach in which excited matter such as a QGP 
consists of a collection of uncorrelated quarks an gluons 
acting as scattering centers for the high energy partons, 
has difficulties accounting for the large jet losses reported 
at RHIC. This is the primary motivation for considering 
synchrotron-like radiation in this work. 







FIG. 3. The l.h.s. shows three diagrams illustrating 
QCD bremsstrahlung radiation: the small circle represents 
the source of the perturbative field. The r.h.s. shows three 
similar diagrams illustrating QCD bremsstrahlung radiation 
(a,b) and synchrotron- like radiation (c). The strong chromo- 
magnetic field is assumed to be in between two thin vertical 
lines, where the full (dressed) propagators are indicated by 
thicker lines. 



B. QCD bremsstrahlung 

For completeness we start with a brief reminder of or- 
dinary QCD bremsstrahlung. This effect was first con- 
sidered by Bertsch and Gunion |2fJ, who calculated the 
three diagrams shown in the l.h.s of Fig. || for soft gluon 
emission. Adding their squares Q at high collision ener- 
gies relative to both the momentum transfer l t from the 
target parton and also the transverse momentum of the 



"The interference terms were also discussed in that paper. 
They are related to form factors which are relevant to radia- 
tion inside hadrons, which we are not interested in here. 



radiated gluon <f t , they have shown that the number of 
gluons emitted is 



dN n 



;2 



dy<Pq t ^ q2(q t -l t f 



(38) 



where C 2 is a color factor (3 for qq scattering with gluon 
emission). Note that at zero momentum transfer l t — 
the radiation vanishes. Also, when the denominators 
become small, they generate two cones of radiation, along 
the initial and final direction of the jet. Integrating this 
result over the rapidity y and the transverse momentum 
qt of the gluon, yields the standard logarithmic factor. 

In our case the field crossed by a jet parton is classical 
and non-perturbative with A ~ l/g, and the radiation is 
described by the modified diagrams shown in the r.h.s. 
of Fig. ||. The motion inside the field is described by 
fully dressed propagators, and the power radiated from 
the slab has been described above. Assuming the chro- 
momagnetic field to be confined to a finite slab, requires 
that we also add the bremsstralunhg diagrams from the 
in and out motion as illustrated in the first two dia- 
grams of Fig. |. The standard QCD variant of the QED 
Weizsacker- Williams (WW) approximation can be used, 
with the so called "parton-in-parton" DGLAP splitting 
functions [fill 



dlrifi 



2n 



(39) 



where x is the standard parton momentum fraction and 
i,j=q,g. The splitting functions are known for all values 
of x, and the total energy loss for quark/gluon due to 
splitting is fill] 



-1^0.28a s \n(Ql igh /Qi w ) 



« 1.0a. 



In (Q 2 



high 



IQ 



2 \ 

low/ 



(40) 



(41) 



where the second integral was regulated by setting 
a^max = m*/P, taken to be 0.95 in the estimate. The 
parton-to-parton splitting happens twice, i.e. in and out. 
On the way in, the 2 scales that define the DGLAP evo- 
lution of the jet are the kick in the scattering process and 
the parton virtuality while hitting the magnet. On the 
way out, the 2 scales are the kick from the slab and the 
final scale of the parton in matter or its hadronization 
scale (whichever is larger). 

We close this summary section by the following qual- 
itative comments: i. The strong bending of the partons 
happens rarely, so it is not subject to the LPM effect; 
ii. The in and out bremsstrahlung effect depends weakly 
(logarithmically) on the field strength provided it is large 
enough to allow for the separation of the 2 cones of ra- 
diation; iii. The magnitude of the relative energy loss 
by bremsstrahlung alone is of order AE/E ~ 1/4, and 
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maybe alone sufficient to explain the expected energy 
loss for gluon jets at RHIC; iv. The contribution from 
diagram (a) is in general small if the jet has been just 
produced in a hard collision, since the virtual field does 
not have sufficient time to form. 



C. Jet quenching in the Color Glass Condensate 

If the initial excited glue is not a set of incoherent glu- 
ons with occupation numbers n ~ O(l), as it is the case 
in an equilibrated QGP, but a coherent classical field with 
large occupation numbers n ~ 0(\/g 2 ), the radiation 
losses are synchrotron-like. Coherence helps, because all 
coherent quanta work in the same direction in space and 
color space, providing larger acceleration. 

To assess the amount of synchrotron radiation in the 
CGC phase, we recall some useful numbers from the nu- 
merical analysis of the SU(3) version of the CGC carried 
in pfij] . At RHIC energies the saturation scale Q s was 
found to be 1.3 GeV. The initial classical CGC field was 
found to abelienize in a time tcgc of order Q s tcgc ~ 3. 
In this regime, the gluon energy density was found to be 
e/Qg ~ 0.17/g 2 with an approximatly thermal momen- 
tum distribution. The transverse energy per quantum 
was found to be 1.66 Q Sl resulting in an effective tem- 
perature of 1 GeV. This is of course only apparent as 
the underlying evolution is classical and originates from 
a coherent state. The field strength F (the r.m.s. com- 
bination of electric and magnetic fields) is about 



gF~ 0.58 Qf. - 1 GeV 



(42) 



giving a quark with e a = g\/3/6 a kick of order e a Fr ~ 
0.66 GeV, and about twice that for a gluon. 

Before substituting these numbers into our expressions 
for the synchrotron radiation loss, we need to do some 
relevant color sums 



C A 



(43) 



The values for "penetrating gluons" of kind 3 and 8 
are 0.30 and 0.22, respectively. For all gluons, we have 
J2a=i 8 — 2.06, with about 1/2 originating from the 
undeflected gluons and 1/2 originating from the deflected 
ones (of course with the extra penalty factor at low uj as 
explained above). For the estimate to follow we will use 
a color factor of 1.5. 

Using the CGC numbers just quoted, we find that the 
relative energy loss of a quark by synchrotron radiation 
in a time tcgc is 



AE, 



CGC 



E 



0.3 



H 



lGeV" 



2/3 / Atcgc \ 
\0.5fm) 



lGeV\ 1/3 



D. Jet Quenching on the Exploding Sphalerons 

The energy and the Chern-Simons number of the re- 
leased (turning) coherent state is 

E = ^£ 
Aap 



(45) 



In so far, the dimensionless parameters have been deter- 
mined in two ways. First, by minimizing the QCD poten- 
tial for fixed Chern-Simons number with the result fllOll 



,2\2 



£ = (1- K\ 

AA = sign( K )(l-| K |) 2 (2+|K|)/2 



(46) 



Eliminating n yields the potential profile E(Nqs)- Sec- 
ond, by maximizing the partial parton-parton cross sec- 
tion with the result [fl2| 



£ = (E/M s ) 
N = (E/M s ) 2 ^ 5 



(47) 



At the sphaleron point the partial cross section is max- 
imum, with £ = J\f = 1. Only this case will be con- 
sidered here. Using instanton vacuum physics we ob- 
tain a sphaleron mass Ms = 3 GeV and a sphaleron size 
p = 1/3 fm. 

Upon release in Minkowski space, the sphaleron state 
evolves classically in real time through the classical Yang- 
Mills equations, as was originally done in the electroweak 
theory (27j. For SU(2) Yang-Mills, this evolution was 
recently carried out in JlOfl both numerically and analyti- 
cally. The turning states were found to explode into thin 
shells of coherent gluonic fields. For our purposes, we 
just recall that the shell for t,r 3> p has a very simple 
radial energy density 



47rr 2 e(r, t) 



8tt 

2 2 

9 P 



p 2 + (r 



(48) 



At large time t ^S> p the corresponding gauge field is 
purely transverse, with equal chromoelectric and chromo- 
magnetic fields. The prompt sphaleron configuration re- 
leased in parton-parton scattering carries initially a very 
strong chromomagnetic field, 



H 



f2M & 



\ P 3 



1/4 



1 GeV 



(49) 



E 



The gluon loss is about twice the quark loss. 



In the early phase of the prompt process in heavy ion 
collisions, the escaping sphalerons form a dilute gas. So 
unlike the CGC they cannot affect most of the jets ini- 
(44)tially for times t ~ p. They do affect them as they ex- 
pand into exploding shells. The net synchrotron radia- 
tion loss involves also the transverse density of sphalerons 
per unit rapidity ns and their typical collision volume 
cr(£) dt. Specifically, 
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AE = P(t)n s a(t)dt 



dtt-^ 3+2 



(50) 



where we used the cross section a ~ t 2 and the radia- 
tion loss P ~ H 2 / 3, ~ t~ 2 / 3 . The result formally diverges 
for large times. However, the above reasoning is only 
valid till the single shell expansion remains coherent. As 
we now show, the originally dilute gas of shells quickly 
evolves into a foam-like structure for times t ~ 2-3p pro- 
viding a natural cutoff in the time integral. 

Recent estimates of the number of clusters produced in 
pp and AA collisions as a function of the collision energy 
and centrality, are still rather uncertain. The theoretical 
calculations of the cross section such as [[Til are carried to 
only exponential accuracy, while phcnomcnological stud- 
ies such as [ p8[ have only resulted into an upper estimate. 
Assuming the whole growth with y/s of the pp cross sec- 
tion to be due to the release of 1 sphaleron, and ignoring 
nuclear shadowing, about S = 400 sphalerons are pro- 
duced in central AuAu collision at y/s — 136 GeV at 
RHIC. Since the total rapidity interval is Ay = 4, this 
amounts to about 100 sphalerons per unit rapidity. The 
transverse density of sphalerons per unit rapidity in cen- 
tral AuAu collisions is 



ns 



nR 2 Ay 



1 f m 



(51) 



which results into a foam-like structure for the exploding 
sphalerons for times t ~ 2p ~ 0.7 fm. This time increases 
as the cubic root of the transverse density for smaller 
densities. For a jet piercing a wall with the energy density 
( ft8| ) the integrated kick is about 0.5 GeV for a quark with 
e a = gVi/ft, and about 1 GeV for a gluon. 

Substituting the expression for the field strength of the 
exploding shell into the synchrotron radiation loss 
(|3l|), yields J F 2 ^ 3 dt — F^/J^irp. For two overlapping 
shells we add the fields in quadrature and estimate that 
the maximal field -F max ~ 0.2 GeV 2 in the foam phase 
at time ifoam ~ 2p. Substituting all this into the syn- 
chrotron energy loss formulae, we finally obtain for the 
quark loss 



AE 
~E~ 



0.21 



H \ 2/3 /lGeV\ 1/3 



0.2 GcV^ 



E 



(52) 



The gluon loss scales with the pertinent color Casimir 
and is about twice larger. 



IV. SUMMARY AND DISCUSSION 



angle scattering subject to the Landau-Pomeranchuck- 
Migdal suppression, we have synchrotron-like QCD radi- 
ation. The radiation is enhanced by the coherent clas- 
sical fields, providing larger acceleration and radiation 
compared to independent quanta. 

There are currently two mechanisms for the formation 
of strong and prompt classical color fields in heavy ion 
collisions. First, the color glass condensate (CGC) Q] 
is a classical Weiszacker- Williams field of virtual gluons 
initially part of the wavefunction of the colliding nuclei, 
that is made real by the collision. Second, the exploding 
sphaleron-like clusters which are the remnants of (singu- 
lar) instantons in the QCD vacuum. The clusters are not 
only coherent, but evolve into thin shells of strongly local- 
ized fields and become foam-like. Any jet has a pobability 
of about 1 to interact with walls of exploding clusters. 

We have shown that synchrotron radiation loss from 
a 6 GeV quark jet is about 0.17 in the CGC and about 
0.12 in the exploding sphalerons. Gluon jet losses are 
about twice larger. The corresponding bremsstrahlung 
radiation loss on the in and out motion results also in 
a loss of a similar magnitude. At larger jet energies, 
bremsstrahlung loss dominates over synchrotron loss, al- 
though the latter is found to decrease slowly with increas- 
ing energy (as l/E 1 / 3 ). All in all, the mechanisms of jet 
quenching considered here provide in total about 20-30% 
loss, which is about consistent with the empirically re- 
ported jet quenching at RHIC. 

The present considerations of jet quenching in the 
context of the color glass condensate or the exploding 
sphalerons is rather schematic, with quantum effects car- 
ried only to leading order. As emphasized above, how- 
ever, what really matters is not the average radiation 
loss for a jet, but rather the probability for a jet to 
escape without losses. To assess this, detailed simula- 
tions with realistic nuclei geometry are needed. Also, as 
suggested in 11 1, it is possible to experimentally mea- 
sure whether the radiated gluons are produced in a nar- 
row cone around the jet (bremmstrahlung) or not (syn- 
chrotron). One can also do tagged jets by measuring 
photon-jet correlations, and search for acomplementarity 
and pt disbalance. Clearly much more theoretical work is 
needed to make the present estimates more quantitative. 
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The chief idea of this work is that if the initial stage 
of a heavy ion collision produces a coherent classical field 
rather than an a gas of incoherent quanta, one should re- 
consider the theory of all prompt processes (Drell-Yann, 
photons, dileptons, heavy quarks, ...), including the cur- 
rent theory of jet energy losses. Instead of multiple small 



APPENDIX A: CLASSICAL REDUCTION 

In this Appendix we give some helpful steps leading 
to the classical formulae (^0|) . We will show how to re- 
duce the gluon part, since the quark part follows from 
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Schwinger's paper |ll| to which we refer to. We first ex- 
ponentiate the scalar quark propagator in (19), combine 
it with the already exponentiated gluon propagator and 
change the proper-time variables (Feynman parametriza- 
tion) to obtain 



da 



OO rl 

sds I — 

o Jo cos(g A Hs(l-u)) 



g— ism u g-isHt 



xe 



-is(l-u) k^({tan(g A Hs(l-u))/(g A «))—!) 



for the propagators only. The proper-time hamiltonian 
is 



(k — uT\]j 



(A2) 



with a 'mass-shell' condition dUb/dk ~ leading to fc ~ 
ullfe or fc 2 ~ u 2 to 2 in the classical limit. This saddle 
point relation receives corrections at the quantum level. 

Following jl6| we introduce the key change of variables 
that facilitates the identification with the classical radi- 
ation problem 



2uj 

u=-«l 



(A3) 



where the last inequality will be re laxed through the first 
quantum correction. Using (A3) and the substitution 
k 2 



valid to leading order in u/E (classical), we 
find the classical limit to the gluonic contribution to be 



e -i(ELO A ) 2 T 3 /(2Auj) 

cos (Eloat/(2uj) 



(A4) 



which is the part quoted in (po|). The remaining quark 
part follows exactly Schwinger's argument and will not 
be repeated here. 
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